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I  SUMMARY 

This  technical  report  summarizes  the  results  of  experimental 
and  analytical  investigations  on  multiple  liquid  (water)  Impact 
eroslo^of  1100-0  Aluminum  and  316  ^^ainless  steely  junder  contract 
No.  NASW-1608.  fjji  existing  rotating  disk  facility  was  modified 
so  that  six  specimens  might  be  simultaneously  tested  at  three  dif¬ 
ferent  velocities.  High  speed  motion  pictures  as  well  as  con¬ 
ventional  movie  films  with  stroboscopic  illumination  were  taken 
to  show  the  phenomenon  of  liquid  Jet  Impact  with  the  test  material 
in  a  qualitative  manner. 

Using  this  facility^  the  threshold  Impact  velocity  at  which 
visible  erosion  was  observed  after  ten  million  Impacts  was  de¬ 
termined.  The  threshold  water  hammer  stress  on  the  material  cor¬ 
responding  to  this  threshold  velocity,  was  correlated  with  the  •  . 
high  frequency  endurance  limit  corresponding  to  ten  million  fatigue 
cycles  using  the  magnetostriction  vibratory  apparatus.  The  ratio 
between  the  fatigue  endurance  limit  and  the  threshold  water  ham¬ 
mer  stress  is  I.7  fo^  1100-0  aluminum  and  3*0  for  318  stainless 
steel. 

A  mathematical  theory  of  erosion  has  been  recently  developed 
to  predict  the  rate  of  cavitation  erosion  as  a  function  of  ex¬ 
posure  time.  This  theory  was  extended  to  multiple  liquid  Impact 
erosion.  In  this  theory,  it  is  important  to  know  the  value  of  a 
parameter  called  the  shape  parameter,  a.  In  the  original  theory 


it  was  assumed  that  the  shape  parameter  was  the  same  as  that  ob¬ 
tained  in  fatigue  probability  distribution  (Weibull  distribution) 
However,  the  present  investigations  show  that  the  Weibull  shape 
parameter  as  obtained  from  a  high  frequency  fatigue  test  cor¬ 
responding  to  ten  million  life  cycles  is  about  0.5  whereas  the 
shape  parameter  required  to  predict  the  erosion  rates  is  3.0  for 
aluminum  and  5.0  for  316  stainless  steel.  The  implications  of 
this  result  on  the  erosion  theory  are  fully  discussed.  The 
probable  reasons  for  the  variation  in  the  shape  parameter  are 
also  pointed  out.  The  usefulness  of  the  erosion  theory  in  practi 
cal  applications  is  also  brought  forth. 

An  analysis  of  the  results  obtained  indicates  that  the  peak 
rate  of  erosion  varies  as  the  fifth  power  of  the  velocity  of 
impact  and  the  time  at  which  the  peak  rate  is  observed  varies  as 
the  one-fifth  power  of  the  Impact  velocity. 

Scanning  electron  microphotographs  of  eroded  specimens  both 
by  cavitation  and  by  liquid  Impingement  were  compared.  The  pro¬ 
cess  of  material  removal  seems  to  be  very  similar  in  both  the 
kinds  of  erosion  further  justifying  the  use  of  a  common  theory. 1 
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INTRODUCTION 

The  problem  of  erosion  caused  by  multiple  liquid  Impacts 
can  become  a  serious  one  in  the  development  of  future  space 
nuclear  power  systems.  In  addition  it  is  also  Important  in  wet 
steam  or  vapor  turbines  and  in  rain  erosion  of  aircraft  and 
missiles.  The  objectives  of  the  present  research  program  are 
the  understanding  of  the  mechanism  and  the  quantitative  evalua¬ 
tion  of  two  aspects  of  multiple  liquid  Impact  erosion,  namely: 

(i)  the  initiation  of  erosion  and  (11 )  the  rate  of  erosion. 

During  the  past  several  years  of  intensive  research  in  the  allied 
problem  of  cavitation  erosion,  the  above  two  aspects  were  con¬ 
sidered  in  depth  (References  1  and  2).  Some  of  the  significant 
ideas  generated  through  those  efforts  were  extended  to  the  under¬ 
standing  of  multiple  liquid  Impact  erosion. 

Specifically  the  major  tasks  of  this  research  program  are 
as  follows: 

1.  To  determine  experimentally  the  threshold  velocity 
of  Impact  on  two  materials  namely  1100-0  aluminum  and  316  stain¬ 
less  steel. 

2.  To  determine  the  high  frequency  fatigue  stresses 
for  these  two  materials. 

3.  To  relate  the  water  hammer  stresses  corresponding 
to  the  threshold  velocities  with  the  high  frequency  endurance 
limit  of  these  materials. 
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To  determine  the  rate  of  erosion  as  a  function  of 
exposure  time  for  these  two  materials. 

5.  To  correlate  the  experimental  rates  with  theoreti¬ 
cally  predicted  rates, 

6.  To  evaluate  the  dependence  of  the  rate  of  erosion 
on  the  velocity  of  Impact. 

These  tasks  were  accomplished  by  using  the  jet-impact  erosion 
facility  and  the  magnetostriction  vibratory  apparatus.  The  re¬ 
sults  of  these  Investigations  are  presented  In  this  report. 
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EXPERIMENTAL  FACILITY  AND  TECHNIQUES 


The  Jet  Impact  Erosion  Facility 

The  erosion  test  facility  that  was  already  available  at 
the  beginning  of  this  program  consisted  of  a  rotating  disk  on 
which  two  specimen  holders  were  fixed  at  two  diametrically  oppo¬ 
site  locations.  A  one  half  horsepower  motor  (capable  of  rotating 
at  10,000  rpm  at  no  load)  drove  the  rotating  disk  up  to  a  maximum 
peripheral  speed  of  350  feet  per  second,  causing  the  specimen  to 
impact  two  solid  Jets  l/l6  inch  in  diameter.  The  erosion  was 
caused  by  the  impact  of  the  test  specimens  on  the  solid  Jet  at 
controlled  speeds.  This  facility  was  Improved  further  under  the 
present  program  by  replacing  the  drive  motor  with  a  2|  hp  high 
speed  motor  capable  of  sustained  operation  at  20,000  rpm  and  a 
suitable  variable  transformer  power  supply  for  speed  control. 

A  new  rotating  disk  capable  of  holding  six  test  specimens  at 
increased  speeds  was  designed,  fabricated  and  used  in  the 
present  investigations.  The  maximum  peripheral  speed  achieved 
was  700  fps.  The  number  of  Jets  was  increased  to  six.  Figure  1 
shows  the  facility  as  it  exists  now  after  all  the  Improvements. 
Appendix  A  contains  the  basic  outline  of  the  design  of  the 
rotating  disk. 

Photographing  the  Jet  Impact  Phenomenon 

Two  techniques  were  used  to  photograph  the  Impact  of  the 
Jet  with  the  test  specimen.  The  first  method  employed  the  high 
intensity  strobe  light  to  stop  the  motion  so  that  conventional 
movie  pictures  at  24  frames  per  second  could  be  taken.  The 
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second  method  used  a  Hyeam  high  speed  motion  picture  camera 
capable  of  taking  up  to  8000  frames  per  second.  The  photographs 
provided  qualitative  understanding  of  the  liquid  Impact  phenomenon. 
In  addition^  they  showed  that  jets  were  not  broken  up  before  Im¬ 
pact  and  that  the  spray  after  Impact  did  not  hit  the  neighboring 
specimens.  Typical  photographs  are  shown  In  Figure  2. 

High  Frequency  Fatigue  Testing  Technique 

A  magnetostriction  oscillator  (Figure  3)  was  used  to  produce 
alternating  strains  at  the  node  of  a  resonating  rod.  The  basic 
principle  and  the  practical  aspects  of  these  tests  are  described 
In  References  3  and  4.  The  fatigue  specimen  design  used  In  the 
earlier  studies  had  a  sharp  notch  at  the  node  (Figure  A).  The 
notch  sensitivity  of  the  test  materials  at  high  frequencies  Is  an 
unknown.  In  order  to  eliminate  this  limitation,  a  new  dumb-bell 
shaped  fatigue  specimen  was  designed  and  used  for  all  the  fatigue 
tests  under  this  program. 

The  basic  approach  for  the  design  of  the  fatigue  specimen  Is 
to  use  the  theory  developed  by  Nepplras  (5).  From  this  theory 
one  can  get  the  lengths  of  the  fatigue  specimen  as  shown  by  the 
example  In  Figure  5.  Assuming  an  area  ratio  andix/h,  the  value 
of  Hs./'K  may  be  determined  from  Nepplras  '  theory.  In  order  to 
avoid  the  sharp  corners,  a  circular  arc  fillet  of  radius  R  Is 
used.  The  value  of  R  can  be  calculated  from  simple  geometrical 
considerations.  This  method  of  designing  dumb-bell  shaped 
fatigue  specimens  gave  the  dimensions  within  10  to  I5  percent 
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accuracy.  Then  the  final  adjustments  are  made  by  tuning  experi¬ 
mentally.  The  dimensions  of  a  properly  tuned  fatigue  specimen 
are  shown  In  Figure  6. 

In  all  the  fatigue  studi.es,  theoretical  strain,  given 
by  the  following  formula  (due  to  Nepplras  (5))^  was  assumed  to  be 
the  actual  value; 


e 


A 


CA] 


where 

e  -  the  strain  amplitude  at  the  node, 

I  -  the  displacement  amplitude  at  the  anti-node, 

A  -  the  wavelength  In  the  material,  and 

G  -  Magnification  Factor 

_  Strain  In  stepped  specimen _ _ 

~  Strain  In  uniform  specimen  (without  step). 

The  value  of  G  can  be  calculated  from  Nepplras'  theory.  Generally 
one  would  measure  the  displacement  amplitude  ^  and  the  wavelength 
A,  and  calculate  the  strain  using  a  theoretical  value  of  G.  How¬ 
ever, we  attempted  to  verify  these  calculations  by  measuring  the 
actual  strains  generated  at  the  node  In  spite  of  many  practical 
difficulties.  Fixing  a  suitable  strain  gage  to  withstand  the 
high  frequency  fatigue  strains  became  a  problem.  The  strain  gage 
failed  In  fatigue.  The  adhesive  did  not  hold  properly  under 
fatigue  strain.  Since  the  entire  fatigue  specimen  has  to  be  Im¬ 
mersed  In  a  cooling  bath,  the  strain  gages  had  to  be  potted  with 
water  resistant  coating.  These  problems  were  solved  after  sev¬ 
eral  trials  by  using  different  strain  gage-adhesive  systems.  The 
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successful  one  was  the  type  PAB-03N-12-S9  with  EPY  500  adhesive 
supplied  by  Baldwin-Llma-Hamilton  Electronics,  Inc.  The  entire 
strain  gage  was  potted  with  neoprene  and  vacuum  cured  making  it 
impermeable  to  water. 

The  fatigue  specimens  were  calibrated  statically  by  using 
deadweights  hanging  from  the  specimen.  In  order  to  minimize 
bending,  two  universal  Joints  were  provided  fore  and  aft  of  the 
fatigue  specimens.  Assuming  that  the  static  calibration  holds 
good  for  the  dynamic  measurements  also,  the  dynamic  strains  were 
measured  at  various  displacement  amplitudes.  A  comparison  be¬ 
tween  the  measured  strains  and  the  calculated  strains  Is  shown 
In  Figure  7.  This  effort  proved  that  the  theoretical  formula  Is 
accurate  for  our  design. 
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EXPERIMENTAL  PROCEDURE 


As  discussed  previously,  essentially  the  following  three 
types  of  experiments  were  performed  during  this  investigation: 

1.  Determination  of  Threshold  Velocities, 

2.  Determination  of  the  fatigue  endurance  limit,  and 

3.  Determination  of  Erosion  Rates. 

The  erosion  tests  were  carried  out  In  the  Jet  Impact  erosion 
facility  and  the  fatigue  tests  were  conducted  In  the  magneto¬ 
striction  vibratory  apparatus.  The  Important  details  of  the  ex¬ 
perimental  procedure  followed  are  described  below. 

Determination  of  Threshold  Velocities 

Six  test  specimens  were  attached  at  three  radial  locations 
in  the  rotating  disk  described  previously.  This  enabled  the 
observation  of  two  specimens  at  a  preselected  test  velocity  en¬ 
suring  the  reproducibility  and  reliability  of  the  experimental 
observation.  The  test  specimens  were  3/8  Inch  In  diameter.  The 
specimen  surface  was  finished  to  32  rms  or  better  by  hand  polish¬ 
ing,  thus  eliminating  the  possibility  of  mechanical  work  hardening 
affecting  the  test  results.  Such  carefully  prepared  test  speci¬ 
mens  were  attached  to  the  rotating  disk  and  run  at  preselected 
speeds.  Observation  of  the  specimens  to  determine  when  denting 
or  fracture  occurred  was  adjusted  to  suit  the  test  sequence. 

When  tests  were  conducted  at  higher  velocities,  damage  could  be 
observed  In  a  short  time  and  hence  the  specimens  were  observed  at 
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short  Intervals  (every  few  minutes),  M  lower  velocities, 
damage  would  not  occur  for  many  hours;  consequently  observations 
were  made  every  half  or  one  hour.  The  observation  consisted  of 
removing  a  test  specimen  and  of  observing  the  surface  exposed  to 
liquid  impact  with  a  10  X  magnifying  glass  under  side  lighting 
which  would  make  the  initial  indentations  "stand  out"  when  il¬ 
luminated  at  the  proper  angle  to  the  surface.  The  time  taken 
for  the  initiation  of  permanent  plastic  indentations  on  the  sur¬ 
face  of  the  test  specimen  through  this  procedure  was  recorded  at 
different  test  velocities.  The  number  of  impacts  were  calculated 
from  the  number  of  revolutions  made  during  that  time  multiplied 
by  two  for  the  two  impacts  per  revolution.  The  relationship  be¬ 
tween  the  number  of  impacts  and  the  velocity  was  plotted  and  the 
velocity  corresponding  to  10  million  impacts  was  obtained  from 
such  a  plot.  This  velocity  is  designated  as  the  threshold  ve¬ 
locity. 

After  the  appearance  of  plastic  dents  on  the  specimen  sur¬ 
face,  this  area  of  denting  will  continue  to  work  harden  as  the 
test  progresses  resulting  in  the  fracture  of  small  fragments  of 
material  from  the  surface.  In  all  our  experiments,  the  criterion 
for  threshold  was  the  appearance  of  detectable  indentation  with 
the  help  of  a  10  X  magnifier  under  side  lighting.  The  appearance 
of  the  surface  after  10  million  impacts  at  I50  fps  on  316  stain¬ 
less  steel  is  shown  in  Figure  8.  The  surface  roughness  profile 
after  10  million  impacts  at  threshold  velocities  is  shown  in 
Figure  9. 
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Determination  of  High  Frequency  Endurance  Limit 


The  fatigue  specimens  were  vibrated  at  l4.2  kcs  at  controlled 
amplitudes  In  the  magnetostriction  vibratory  apparatus.  The  speci¬ 
mens  were  cooled  by  immersion  in  a  constant  temperature  water  bath. 
The  amplitude  was  obtained  from  the  voltage  output  df  a  precali¬ 
brated  voice  coll.  The  strain  was  calculated  from  the  amplitude 
(using  equation  [A])  as  discussed  previously.  The  stress  was 
obtained  by  multiplying  the  strain  with  the  modulus  of  elasticity. 
The  modulus  of  elasticity  for  the  test  material  can  also  be  de¬ 
termined  with. the  vibrator  by  measuring  the  wavelength  and  the 
frequency  which  will  give  the  speed  of  sound.  The  frequency  was 
accurately  measured  with  a  frequency  counter.  The  time  to  fail¬ 
ure  was  noted  and  the  number  of  cycles  to  failure  was  computed 
from  the  frequency. 

Determination  of  Rate  of  Erosion  as  a  Function  of  Exposure  Time 

The  rate  of  weight  loss  was  determined  periodically  by  weigh¬ 
ing  the  same  specimen  after  exposure  to  impact  at  preselected 
velocities.  Care  was  taken  to  adopt  a  uniform  procedure  in  clean¬ 
ing  and  drying  the  specimen  before  and  after  each  test.  The  rate 
of  volume  loss  was  obtained  by  dividing  the  rate  of  weight  loss 
by  the  density  of  the  material.  The  mean  depth  of  erosion  was 
calculated  by  dividing  the  volume  of  erosion  by  the  area  of 
erosion.  The  eroded  area  was  taken  to  be  equal  to  the  impact 
area  of  the  Jet.  The  actual  eroded  area  can  also  be  measured 
accurately  with  a  plahlmeter. 
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RESULTS  AND  ANALYSIS 


Correlation  of  Threshold  Water  Hammer  Stresses  with  High 
Frequency  Endurance  Limit 

Figure  10  shows  the  relationship  between  the  velocity  of 
impact  and  the  number  of  Impacts  after  which  the  Indentations 
were  observed  both  on  1100-0  aluminum  and  on  316  stainless  steel. 
Several  repeated  observations  showed  that  the  threshold  veloci¬ 
ties  corresponding  to  10  million  impacts  were  50  fps  for  1100-0 
aluminum  and  I50  fps  for  316  stainless  steel  . 

When  a  cylindrical  column  of  liquid  impinges  on  the  surface 
of  a  material,  the  maximum  stress  developed  at  the  point  of  con¬ 
tact  by  the  Impact  (generally  known  as  "water  hammer"  stress)  is 
derived  by  de  Haller  (6)  as 


where 


a,  = 


P  C 
■^m  m 


1  + 


Uj  =  Impact  velocity, 

=  Density  of  Liquid, 

=  Density  of  Material, 

=  Velocity  of  Sound  in  Liquid,  and 
C^  =  Velocity  of  Sound  in  the  Material. 


[1] 
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For  the  present  Investigations  In  which  water  was  the  test  liquid 

and  common  metals  are  the  test  materials,  the  ratio 

Is  small  compared  to  unity.  Then  the  water  hammer  stress  becomes 

<’1  =  ‘21 

If  we  assume  that  the  Impact  stresses  must  exceed  the  fatigue 
endurance  strength  of  the  material  In  order  for  one  to  see  any 
permanent  deformation,  then  the  ratio  of  these  stresses  at  thresh¬ 
old  must  equal  unity.  However,  an  analysis  (Reference  l)  of  some 
Isolated  experimental  results  of  several  earlier  Investigators 
showed  that  the  ratio 


131 

One  of  the  major  alms  of  the  present  Investigations  Is  to 
check  the  a  cove  result  by  experimentally  determining  the  thresh¬ 
old  Impact  stress  and  the  endurance  limit  for  two  materials. 

The  values  of  and  for  water,  were  obtained  from 
published  literature  (for  example  Reference  7)  ^nd  the  values  of 
were  obtained  from  Figure  10.  Hence  the  Impact  stress  P^C^U^ 
could  be  calculated. 

The  high  frequency  fatigue  strength  of  the  two  test  ma¬ 
terials  were  determined  as  described  previously.  The  relationship 
between  the  stress  amplitude  and  the  number  of  cycles  to  failure 
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Is  shown  In  Figures  11  and  12.  The  endurance  limit  at  10  mil¬ 
lion  cycles  is  5000  psl  and  30,000  psi  for  1100-0  aluminum  and 
316  stainless  steel  respectively. 

The  pertinent  data  required  for  analysis  is  shown  in  Table  1. 
It  can  be  seen  that  the  ratio  between  the  high  frequency  endurance 
limit  and  the  impact  stress  is  in  the  vicinity  of  two.  Actually 
it  is  1.7  for  1100-0  aluminum  and  3-0  for  316  stainless  steel. 

This  variation  may  be  attributed  to  several  factors  such  as  the 
following: 

1.  The  detection  of  Initiation  of  damage  is  by  visual 
observation.  This  may  introduce  some  uncertainties  in  the  sur¬ 
face  appearance  of  the  individual  metals  studied. 

2.  Near  threshold  velocities,  the  effect  of  corrosive 
influences  may  produce  surface  discoloration  and  give  rise  to 
differences  in  observed  threshold  velocities.  However,  the  sur¬ 
face  roughness  profiles  rule  out  this  possibility. 

3.  It  is  quite  possible  that  some  other  mechanical 
property  is  responsible  for  producing  the  variation. 

4.  The  de  Haller  equation  to  predict  the  water  hammer 
pressure  on  the  material  may  not  be  accurate  enough  for  the 
present  experimental  conditions,  since  it  is  based  on  one-dimen¬ 
sional  analysis  for  the  impact  of  two  solid  bodies. 
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Effect  of  Time  of  Exposure  to  Multiple  Liquid  Impacts  on 
The  Rate  of  Erosion 

The  importance  of  exposure  time  on  the  rate  of  erosion  was 
recognized  and  highlighted  in  a  series  of  recent  publications 
(8^9^10^11  and  12)  based  mainly  on  cavitation  erosion  research. 
The  typical  relationship  between  the  test  duration  and  the  rate 
of  erosion  in  a  cavitation  erosion  test  is  shown  in  Figure  13 . 
This  relationship  was  divided  into  four  zones.  Figure  l4  as 
follows: 

1.  Incubation  Zone:  The  rate  of  erosion  is  very  small 
during  the  early  part  of  an  erosion  test.  During  this  time,  the 
material  undergoes  permanent  changes  due  to  the  repeated  erosive 
forces.  Hence  •it  is  believed  that  the  erosion  Incubates  during 
this  period. 

2.  Accumulation  Zone:  After  the  incubation,  the 
material  starts  fracturing  due  to  the  repeated  Impacts.  The  rate 
of  loss  of  material  starts  increasing  with  further  exposure  to 
erosion.  Since  this  process  is  similar  to  cumulative  fatigue 
fracture,  it  was  called  energy  accumulation  zone. 

3.  Attenuation  Zone:  Where  enough  material  has  been 
fractured,  the  eroded  surface  becomes  rough  with  visible  deep 
craters.  In  cavitation  erosion  and  in  liquid  impact  erosion, 
the  liquid  covering  these  craters  cushions  the  Impact  pressures 
transmitted  to  the  material.  The  decrease  in  erosion  rate  is 
believed  to  be  due  to  this  attenuation  process  and  hence  this 
zone  is  called  the  attenuation  zone. 
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4.  Steady  State  Zone;  It  Is  experinientally  observed 
that  after  the  attenuation  zone,  the  rate  of  erosion  becomes  very 
nearly  independent  of  the  exposure  time  and  hence  it  is  called 
the  steady  state  zone. 

Although  the  above  results  on  the  effect  of  exposure  time 
were  observed  only  for  cavitation  erosion^  there  were  some  in¬ 
dications  in  the  literature  that  these  effects  were  also  detectable 
in  liquid  impact  erosion  tests  such  as  steam  turbine  erosion  (13^ 
l4)j  Jet  impact  erosion  (15 )  and  rain  erosion  (l6).  A  summary  and 
discussion  of  all  these  results  with  a  mathematical  analysis  was 
presented  by  Heymann  (17). 

At  this  Juncture^  it  became  crucial  to  conduct  a  few  sys¬ 
tematic  experiments  to  verify  whether  such  effects  were  noticeable 
in  the  case  of  multiple  liquid  Impact  erosion.  As  shown  in  Fig¬ 
ures  15  through  18,  these  effects  are  Indeed  present.  The  ma¬ 
terials  used  for  the  test  were  316  stainless  steel  and  1100-0 
aluminum.  The  results  shown  correspond  to  six  velocities  for  each 
material.  It  may  be  noted  that  the  incubation  period  is  very 
noticeable  at  lower  speeds  whereas  the  peak  rate  of  erosion  is 
very  pronounced  at  higher  Impact  speeds. 

It  is  necessary  to  understand  and  predict  these  non-linear 
time  effects  quantitatively  in  order  to  achieve  meaningful  cor¬ 
relations  in  the  laboratory  and  to  extrapolate  laboratory  data  to 
field  systems.  With  this  aim  in  mind,  a  mathematical  theory  of 
erosion  has  been  developed  recently  (2)  and  is  summarized  briefly 
in  Appendix  B. 


Correlation  with  the  Theory  of  Erosion 

As  shown  In  Appendix  B,  the  following  normalized  differential 
equation  of  erosion  can  be  derived  with  the  two  assumptions  In¬ 
volving  "accumulation"  of  energy  and  the  "attenuation"  of  energy 
and  with  the  definition  of  Intensity  of  erosion: 


2n+l 


dl  .  r  I 

"T —  +  k  r~7 
dr  -  1/n 

T) 


1^.0 
-  dT 


where 


I  =  I  /I  , 

e'  max 

1  =  Maximum  Intensity  of  erosion, 

max 

T  =  t/ti, 

I  =  Intensity  of  erosion  at  any  time  t, 
e 

ti  =  the  time  corresponding  to  the  maximum  Intensity 
of  erosion, 

•Q  =  the  efficiency  of  erosion  at  any  time  t, 

r|i  =  bhe  efficiency  corresponding  to  the  peak  In¬ 

tensity  of  erosion, 

k  -  ,  and 

n  =  the  attenuation  exponent. 
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The  general 
tlon  Is  found  to 


solution  of  this  normalized  differential  equa¬ 
te 


f 


1  +  k2±i 

n 


TldT 


n/n+1 


[5] 


In  order  to  calculate  I  as  a  function  of  the  value  of  n 
and  the  function  Ti(t)  must  be  known.  Prom  the  experience  of 
underwater  explosions^  It  Is  known  that  the  shock  pressure  at¬ 
tenuates  Inversely  as  the  distance  travelled  (I8).  Since  the 
Intensity  varies  as  the  square  of  the  shock  pressure,  the  value 
of  n  may  be  assumed  to  be  2. 

In  the  original  formulation  of  the  elementary  theory  (2), 
the  nature  of  the  function  ri(t)  was  assumed  to  be  similar  to  the 
fatigue  probability  function.  The  physical  basis  for  this  as¬ 
sumption  may  be  enunciated  as  follows:  A  fractured  particle  has 
absorbed  energy  from  the  Impact  forces  over  a  number  of  cycles. 
After  the  cumulative  absorption  of  these  Impact  forces,  the 
particle  fractures  from  the  specimen  thereby  producing  a  loss 
of  material.  In  this  sense,  the  efficiency  of  the  erosion  pro¬ 
cess  Is  associated  with  the  fracture  of  a  particle  of  the  ma¬ 
terial.  If  we  consider  the  particle  as  a  fatigue  specimen,  then 
the  probability  of  failure  of  the  particle  after  sometime  Is  a 
statistical  function.  Based  on  this  approach.  It  may  be  Inferred 
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that  the  efficiency  of  energy  absorption  (the  function  Ti(t))  Is 
associated  with  the  probability  of  failure  of  a  particle  after 
a  given  time.  There  are  several  statistical  distribution  func¬ 
tions  advanced  in  the  literature  to  represent  the  probability  of 
fatigue  failure.  Specifically,  the  Weibull  distribution  (19) 
was  chosen  for  the  analysis  because  of  Its  wide  applicability. 
Hence  a  Weibull  type  distribution  function  of  the  following  form 
was  assigned  to  t)  : 

T)  =  1  -  exp(-T*^)  [6] 

The  parameter  a  Is  called  the  Weibull  shape  parameter  and 
it  depends  upon  the  material  as  well  as  the  stress  level.  As 
can  be  seen,  the  above  distribution  becomes  a  simple  exponential 
distribution  where  a  =  1  and  the  Rayleigh  distribution  when 
a  =  2.  It  corresponds  to  a  nearly  normal  distribution  when 
a  =  3.57  (20).  Using  different  values  of  a  In  Equation  [6],  the 
solutions  to  Equation  [5]  were  evaluated  and  shown  In  Figure  19 
as  a  plot  of  the  relative  Intensity  against  the  relative  ex¬ 
posure  time  with  a  as  the  third  parameter. 

At  this  juncture,  the  following  points  about  the  elementary 
theory  discussed  above  are  worthy  of  consideration; 

1.  The  theory  is  not  llm.lted  to  the  use  of  Weibull 
type  distributions  to  represent  the  efficiency  of  energy  absorp¬ 
tion.  In  fact,  it  is  also  possible  to  use  a  direct  relationship 
between  the  energy  absorption  capacity  and  the  number  of  cycles 
of  impacts.  As  of  now  no  such  quantitative  relationship  exists 


19 


although  one  may  develop  an  empirical  relationship  to  fit  the 
available  experimental  data.  For  example  Manson  (21)  has  pre¬ 
sented  experimental  evidence  to  show  how  work  hardening  or  work 
softening  of  metals  proceeds  with  number  of  fatigue  cycles.  The 
fact  that  the  phenomenon  of  work  hardening  or  work  softening  Is 
Important  in  the  erosion  process  has  been  pointed  out  by  many 
investigators  including  Engel  (22).  However,  the  state  of  the 
art  is  not  adequate  to  develop  a  general  mathematical  relation¬ 
ship  to  fit  the  experimental  data.  Because  of  this  reason,  the 
Weibull  type  distributions  were  used  as  a  first  attempt. 

2.  In  the  elementary  theory  only  the  efficiency  func¬ 
tion  Ti(t)  was  assumed  to  be  a  statistical  distribution.  However, 
the  drop  size  distribution  (or  bubble  size  distribution  In  the 
case  of  cavitation)  and  the  surface  roughness  distribution  are 
both  statistical  quantities.  For  this  reason,  it  was  pointed  out 
in  Reference  2  that  a  general  theory  of  erosion  would  lead  to  an 
erosion  distribution  function,  which  wcSuld  take  into  account  all 
of  the  above  parameters. 

3-  Again,  It  was  assumed  in  the  elementary  theory, 
that  the  value  of  the  shape  parameter  a,  was  exactly  the  same 
as  one  would  obtain  from,  a  fatigue  life  distribution  for  the 
material.  Another  assum.ption  Involved  Is  that  the  value  of  a 
is  unique  at  all  stress  levels.  However,  It  is  known  that  the 
value  of  the  shape  parameter  depends  upon  the  stress  level,  the 
presence  of  notches,  the  environment  and  the  material  (19). 
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With  this  background j  let  us  analyze  the  experimental  data 
shown  in  Figures  15  through  l8.  These  data  are  reduced  in  the 
non-dimensional  form  by  dividing  the  rate  of  loss  at  any  time  by 
the  peak  rate  of  loss^,  the  latter  being  obtained  from  the  curve 
which  fitted  the  experimental  data  best.  The  time  is  normalized 
with  respect  to  the  time  at  which  the  peak  erosion  rate  is  ob¬ 
served.  The  data  points  collected  at  five  different  velocities 
are  shown  in  Figures  20  and  21  for  1100-0  aluminum  and  for 
316  stainless  steel  respectively.  The  solid  lines  in  these 
figures  are  theoretical  predictions  of  Equation  C5]*  Weibull 
distributions  of  the  type  given  by  Equation  [6]  are  used  for  the 
efficiency  function  assuming  a  =  3  for  1100-0  aluminum  and  a  =  5 
for  316  stainless  steel.  This  result  in  itself  is  very  useful  in 
predicting  the  erosion  rates.  However^  it  also  leads  to  the 
following  question  of  scientific  interest. 

How  does  the  value  of  the  shape  parameter,  a  as  determined 
in  the  above  erosion  test  compare  with  the  shape  parameter,  a  as 
obtained  from  an  actual  fatigue  life  distribution  for  each 
material?  In  order  to  make  a  distinction  between  the  two  shape 
parameters,  let  us  call  the  shape  parameter  for  the  erosion 
distribution  as  hi,  and  the  shape  parameter  for  the  fatigue  life 
distribution  as  hg . 

The  fatigue  life  distribution  for  the  two  materials  were 
experimentally  determined  using  the  high  frequency  fatigue  tech¬ 
nique  described  earlier  in  this  report.  Twenty-nine  samples  of 
1100-0  aluminum  were  tested  at  a  stress  level  of  5OOO  psi  cor¬ 
responding  to  a  mean  life  of  10  million  cycles  (see  Table  2). 


Also  shown  in  this  table  are  the  Welbull  distribution  func¬ 
tions.  The  data  reduction  has  been  carried  out  as  suggested  in 
Reference  2.  Figure  22  shows  the  plot  of  the  results  in  a  Welbull 
probability  paper.  The  Welbull  shape  parameter,  actually 

varies  with  the  number  of  cycles  particularly  below  500,000  cycles. 
However,  for  the  most  part,  the  Welbull  shape  parameter,  012,  is 
nearly  0.5.  The  mean  life  is  3  x  lo'^  cycles.  Similar  results  are 
shown  in  Table  3  and  in  Figure  23  for  316  stainless  steel.  The 
value  of  Qg  is  close  to  O.A  for  this  case.  Again,  the  Welbull  plot 
does  not  seem  to  be  a  straightline  for  the  entire  range  of  the 
life  cycles.  For  example  in  the  case  of  316  stainless  steel  data, 
at  the  lower  end  of  the  life  cycles,  the  value  of  cfe  is  close  to 
1.0.  The  same  trend  is  noticeable  for  the  case  of  1100-0  aluminum. 

In  any  case,  it  is  clear  that  the  value  of  the  shape  param¬ 
eter  as  as  obtained  by  high  frequency  fatigue  tests  on  dumb-bell 
shaped  fatigue  specimens  at  a  stress  level  corresponding  to  the 
mean  life  of  twenty  to  thirty  million  cycles  is  much  lower  than 
the  value  of  the  shape  parameter  ai,  as  determined  from  the 
erosion  curve.  The  following  factors  must  be  taken  into  con¬ 
sideration  in  understanding  the  above  result. 

1.  The  erosion  was  produced  at  much  higher  stress 
levels  than  the  stress  level  at  which  fatigue  distribution  was 
determined.  It  is  known  that  the  value  of  a  is  much  higher  at 
higher  stress  levels  (Reference  I9). 
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2.  The  fracture  in  erosion  is  produced  at  sharp 
corners  and  stress  concentrations  whereas  the  dumb-bell  shaped 
fatigue  specimen  is  rather  idealized  when  compared  to  the  fractured 
particle. 

3.  As  pointed  out  previously,  the  influence  of  other 
factors  such  as  surface  roughness  distribution  and  drop  size 
distribution  must  also  be  considered  in  addition  to  the  fatigue 
life  distribution. 

While  these  aspects  are  under  continued  investigations,  the 
results  so  far  obtained  do  provide  very  useful  practical  results. 
For  example,  the  effect  of  test  duration  may  be  quantitatively 
reduced  in  a  generalized  form.  It  can  be  mathematically  expressed 
if  the  shape  parameter  ai  is  determined  through  an  erosion  test 
for  each  material.  As  shown  in  Appendix  B  this  relationship  can 
be  Integrated  so  that  the  cumulative  depth  of  erosion  may  be 
predicted  at  different  Intensities  of  erosion. 

The  Relationship  Between  the  Impact  Velocity  and  the 
Rate  of  Erosion 

The  relationship  between  the  velocity  of  Impact  and  the  rate 
of  erosion  is  very  Important  in  understanding  the  phenomenon  of 
erosion  as  well  as  in  extrapolating  laboratory  data  to  practical 
problems.  Heyraann  (23)  has  summarized  most  of  the  available  in¬ 
formation  on  this  aspect.  Without  considering  the  effect  of  test 
duration,  several  investigators  have  found  that  the  rate  of  ero¬ 
sion  varies  as  some  power  of  the  Impact  velocity  as  listed  below 

(23). 
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Author 


Exponent 


Honegger  2.0 

Pearson  2.6 

Hoff  et.al.  5-1 

Hobbs  5 


After  analyzing  all  these  data,  Heymann  (23)  came  to  the  con¬ 
clusion  that  the  exponent  will  generally  be  of  the  order  of  5. 
Heymann  also  pointed  out  that  the  exponent  is  close  2.5  if  the 
threshold  velocity  is  subtracted  from  the  actual  Impact  velocity. 

In  order  to  see  how  far  these  conclusions  agree  with  our 
experimental  results,  the  peak  rate  of  volume  loss  as  shown  in 
Figures  15  through  I8  are  plotted  as  a  function  of  Impact  velocity 
in  Figure  24,  (See  also  Table  4).  The  solid  lines  in  Figure  24 
correspond  to  a  fifth  power  variation.  In  other  words,  the  peak 
rate  of  erosion  does  vary  as  the  fifth  power  of  the  velocity  of 
Impact,  Similarly  the  time  at  which  the  peak  rate  is  observed 
is  also  shown  plotted  against  Impact  velocity  in  Figure  25  and 
the  solid  lines  in  this  figure  correspond  to  a  one  fifth  power 
variation.  To  summarize  these  results. 
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where 


I  is  the  maximum  Intensity  of  erosion, 

peak 

ti  is  the  time  at  which  the  peak  intensity  is 

observed,  and 

U  is  the  impact  velocity. 

In  Figure  26,  the  peak  rate  of  erosion  is  plotted  against 
the  difference  between  the  actual  impact  velocity  and  the  thresh¬ 
old  velocity  for  the  test  material.  Here  again  the  solid  lines 
correspond  to  an  exponent  of  five  and  may  be  expressed  empirically 
as 

I  .a  (U-U^)^  .  [9] 

peak  '  T' 

where  U,p  is  the  threshold  velocity  for  the  material.  This  re¬ 
sult  does  not  agree  with  the  conclusion  reached  by  Heymann  (23). 
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COMPARISON  OP  LIQUID  IMPACT  EROSION 
AND  CAVITATION  EROSION  THROUGH 
SCANNING  ELECTRON  MICROSCOPY 

An  attempt  has  been  made  In  this  report  to  make  use  of  the 
theory  of  erosion  developed  on  the  basis  of  ideas,  concepts  and 
experience  accumulated  through  Intensive  research  on  cavitation 
erosion.  Hence  it  is  essential  to  consider  how  far  the  material 
removal  process  in  multiple  liquid  Impact  erosion  is  similar  to 
the  material  removal  in  cavitation  erosion.  The  main  objective 
approach  is  to  establish  a  common  quantitative  frame 
work  for  both  the  kinds  of  erosion. 

In  order  to  compare  the  indentations  and  fracture  produced 
by  the  two  different  erosion  processes,  carefully  prepared  and 
metallurglcally  polished  specimens  (see  the  photomicrographs  in 
Figure  27)  were  subjected  to  both  types  of  erosion  separately. 
Cavitation  erosion  was  produced  in  a  magnetostriction  type  vi¬ 
bratory  apparatus  whereas  the  multiple  liquid  impact  erosion  was 
produced  in  the  facility  described  in  this  report.  Scanning 
electron  micrographs  were  taken  with  the  kind  assistance  of  Mr. 
Robert  Anstead  of  Goddard  Space  Plight  Center  and  reproduced  in 
Figure  28.  Since  the  eroded  surface  is  very  rough,  scanning 
electron  micrography  offers  the  most  suitable  depth  of  focus 
which  cannot  be  obtained  either  by  light  microscopy  or  by  con¬ 
ventional  electron  microscopy.  (Scanning  electron  microscopy  is 
a  relatively  new  technique  in  which  the  electrons  liberated  from 
the  specimen  are  scanned  in  synchronism  with  the  electron-probe 
scanning  the  specimen  (24)). 
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As  seen  In  Figure  2S,  the  similarity  of  the  mechanism  of 
material  removal  is  great.  The  plastic  indentations  overlapping 
each  other  and  producing  raised  craters  are  clearly  seen  in  both 
cases.  Such  overlapping  indentations  cause  the  fracture  of 
particles  of  the  test  material.  These  observations  further  Justify 
the  use  of  the  erosion  theory  developed  for  cavitation  erosion  to 
liquid  Impact  erosion. 
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CONCLUSIONS 


The  following  conclusions  may  be  reached  from  the  experi¬ 
mental  and  analytical  investigations  on  the  multiple  liquid  Impact 
orosion  phenomenon  so  far  carried  out: 

1.  The  ratio  between  the  high  frequency  endurance  limit 
and  the  impact  stress  is  in  the  vicinity  of  two.  Actually  it  is 
1.7  for  1100-0  aluminum  and  3.0  for  316  stainless  steel. 

2.  The  rate  of  erosion  caused  by  multiple  liquid  im¬ 
pacts  is  very  much  dependent  on  the  exposure  time.  The  rate  of 
erosion  is  small  in  the  beginning  (incubation  zone),  increases 
(accumulation  zone),  then  decreases  (attenuation  zone)  and  reaches 
a  relatively  steady  rate  (steady  state  zone).  This  behavior  is 
very  close  to  the  observations  made  on  cavitation  erosion. 

3.  The  tJieory  of  erosion  developed  to  predict  cavita¬ 
tion  erosion  rates,  was  used  to  correlate  with  the  rate  of  ero¬ 
sion  produced  by  multiple  liquid  Impact  erosion.  The  shape 
parameter  as  obtained  in  a  high  frequency  fatigue  curve  at  stress 
levels  corresponding  to  ten  million  cycles  is  much  lower  than  the 
shape  parameter  required  to  predict  the  erosion  rates.  However, 
if  a  proper  shape  parameter  is  assumed  by  inspection  of  the 
normalized  erosion  data,  then  the  correlation  between  the  theory 
and  experiments  is  good.  The  Implications  of  these  findings  and 
the  further  modifications  to  the  theory  as  a  result  of  these 
findings  are  fully  discussed  in  the  report. 
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4.  The  peak  rate  of  erosion  varies  as  the  fifth  power 
of  the  velocity  of  impact.  The  time  at  which  the  peak  rate  of 
erosion  is  observed  varies  as  the  one-fifth  power  of  the  Impact 
velocity. 

5.  Scanning-electron-mlcrophotographs  of  erosion 
caused  both  by  cavitation  and  by  multiple  liquid  Impacts  show 
close  similarity  of  the  process  of  material  removal  by  both  kinds 
of  erosion. 
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APPENDIX  A 


DESIGN  OP  ROTATING  DISK  CAPABLE  OP  HOLDING  SIX  SPECIMENS 

Plgure  29  shows  a  sketch  of  the  rotating  disk  capable  of 
holding  six  specimens.  The  specimens  are  located  at  three  radii 
thereby  yielding  test  results  at  three  velocities  for  every  run. 
Since  there  are  two  specimens  at  each  radius^  the  reproducibility 
of  data  is  verified  at  each  run. 

Power  Requirements :  Since  the  rotating  disk  contains  both  a 
circular  disk  and  two  whirling  arms  extending  from  the  disk,  the 
power  required  to  rotate  the  disk  was  calculated  in  two  separate 
steps . 

1.  Power  consumed  by  the  disk  is  given  by 


where 


C  cl'^R  ' 
m2  D 


[Al] 


C^  is  the  moment  coefficients  given  in  Reference  25, 

p  is  the  density  of  air, 

0)  is  the  angular  velocity  of  disk,  and 
Rjj  is  the  radius  of  disk. 

2.  Power  Consumed  by  the  Whirling  Arms  is  given  by 


P 

a 


Ct^p  cd®C 
D^a 


[A2] 
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where 


is  the  drag  coefficient  of  the  aerofoil, 

C  Is  the  chord  length  of  the  aerofoil, 

R  Is  the  outer  radius  of  the  arm,  and 

o 

R.  Is  the  Inner  radius  of  the  arm. 

1 

3.  Power  consumed  by  the  specimens  Is  given  by 

where 

Is  the  drag  coefficient  for  specimen, 

Jjo 

A  Is  the  projected  area  of  the  specimen,  and 
U  Is  the  peripheral  velocity  of  the  specimen. 

Stress  Calculations:  The  area  of  cross-section  of  a  whirling 
arm  In  which  the  stresses  developed  due  to  centrifugal  forces 
are  constant  from  root  to  tip  Is  given  by  Reference  26  as  follows 


where 


exp 


2ga  j 


[A4] 


a  Is  the  root  area, 
r 

a  Is  the  tip  area, 

W  Is  the  weight  per  unit  volume  of  the  material, 
o  Is  the  allowable  stress,  and 
g  Is  the  acceleration  due  to  gravity. 
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similarly  the  following  equation  gives  the  variation  of  thick¬ 
ness  along  the  radius  for  a  rotating  disk  of  constant  stress 
(Reference  26) 


where 


t 

o 

t 


exp 


Wao^r^ 

2ga 


t  Is  the  thickness  at  the  tip, and 
o 

t  Is  the  thickness  at  any  radius,,  r. 


[A5] 


In  addition  to  these  stresses,  the  centrifugal  stresses 
caused  by  the  weight  of  the  specimen  should  also  be  Included. 
Calculations  using  the  above  relationships  proved  thai;  the  con¬ 
figuration  shown  In  Figure  29  could  be  successfully  built  with 
7075  aluminum  alloy  and  operated  up  to  750  fps  with  a  2-1/2  hp 
motor.  Subsequent  operational  experience  shows  that  the  above 
design  is  sound. 
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APPENDIX  B 


DERIVATION  OF  ELEMENTARY  THEORY  OF  EROSION 


Definition  of  the  Intensity  of  Erosion 

The  Intensity  of  erosion  Is  defined  as  the  power  absorbed 
by  unit  eroded  area  of  the  material  (1,2  and  8)  and  Is  given  by 


I 

e 


rS 
_ e 

t 


[Bl] 


where 

r/t  Is  the  rate  of  depth  of  erosion,  and 

S  Is  the  erosion  strength, 
e 

If  the  rate  of  depth  of  erosion  Is  expressed  In  differential 
form  then 


I 


e 


=  S 


e  dt 


[B2] 


where  r  Is  the  mean  depth  of  erosion  as  measured  from  the 
original  surface  of  the  material  (Figure  30). 

Basic  Assumptions 

The  two  basic  assumptions  of  the  theory  are  as  follows: 

1.  The  attenuation  of  the  intensity  of  collapse  or 
collision  is  assumed  to  be  Inversely  proportional  to  the  nth 
power  of  the  radial  distance.  This  may  be  mathematically  stated 
as 
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I. 

1 


where 


(r«c)" 


tB3l 


Is  the  Intensity  of  Impact, 

I  Is  intensity  of  collapse  or  collision, 

L/ 

r  is  the  radial  distance  (Figure  30), 

A  is  a  proportionality  constant  with  length  as 
dimension,  and 

n  is  the  attenuation  exponent. 


2.  The  second  assumption  is  that  the  intensity  of 
erosion  is  related  to  the  intensity  of  impact  in  the  following 
manner 


I 

e 


[B4] 


where  ti  is  a  material  property  governing  the  efficiency  of 
energy  absorption.  It  is  also  further  assumed  that  the  efficiency 
of  absorption,  rj  is  time  dependent 

n  =  Ti(t) 

Derivation  of  the  Differential  Equation  of  Erosion 
Differentiating  Equation  [b4]  we  get 


dl 


^  dl. 

-  =  — ^  I .  +.71  — rr 
1  '  dt 


dt  dt 


[B5] 
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Differentiating  Equation  [B3]j  we  get 


dl 


J.  =  (-n) _ ^ 


[b6] 


Furthermore  from  Equation  [B3] 


n 

,  c 

(r+r^)  = 


(r+rc) 


n+1 


Prom  Equation  [B2] 


Hence  Equation  [b6]  becomes 


dl. 
_ 1 

dt 


n 


=  (-n) 


I. 

1 

o 

H 

< 

i  h  1 

dr 

I 

e 

dt 

S 

e 

c 

n+l/n 


(n+l)/n  S 
n^  \  ^  e 


A  I 


-n  a’^I 


(n+l)/n 


A^I 


n+l/n 


[B7] 


[b8] 


[B9] 
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(2n+l  )/n 


kl 

e 


_e  ^  ^  0 

■q  dt  ^ 


[BIO] 


Normalized  Differential  Equation 

The  Equation  [BIO]  can  be  normalized  with  respect  to  the 
parameters  corresponding  to  the  maximum  Intensity  oT  damage  ob¬ 
served  In  experiments.  At 


36 


t 

=  ti 

I 

e 

=  I 

max 

dl 

e 

dt 

=  0 

T1 

=  Til 

Then 

II 

;  dx  = 

dt 

ti 

I 

f  e 

II 

IM 

• 

dl 

e 

^  I 

max 

I 

max 

1 

II 

J  dii  = 

dr] 

Til 

Since 

dl 

e 

dt 

=  0  at  t  = 

ti 

Equation  [BIO]  becomes 

2n+l 


n 


kl 


max 


l/n 


max  dn 
Til  dt 


t=ti 
- (n+1 )/n 


k  = 


max 


Til 


n-l/n 


dri 

dt 


1  "ti 


[Bll] 
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General  Solution  of  the  Normalized  Differential  Equation 


Let 


z 


dz 


i 


n+l  \ 

/ 


-2n-l 
,  n 

ILtl  I- 
n 


dl 


Equation  [B12]  becomes 


n 

n+l 


2n+l 


2n+l 

_  _  n 

k  I 


-  dT 
1 


0 


dz  n+l  ^ 

dT  n  -  dT  ^ 

T1 


n+l  k 


This  Is  of  the  form 


f 


faCri) 


General  solution  of  this  equation  Is 
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z  =  exp[-/fi  (ri)dT]  /  [exp  /fi  (rj )dT]f2  (rj )dT 


C  exp[-/fi  (fi)dT]  =  exp[-/fi  (■n)dT][P(Ti)  +  C] 


/fl(Tl)dT  = 


I-  n+1  1  dri  , 

J - j  dT 

•'  n  -  dT 


£n  T] 


=  £n(^) 


n+l/r 


-(n+l)/n 

exp[-/fi  (Ti)dT]  =  exp[-in(Ti)  ] 


-  n+l/n 


=  (n) 


exp[/fi  (Ti)dT]  =  (ri) 


(n+1 )/n 


F(n)  -  /  ^  ^  ■3. 

T) 


Hence  the  general  solution  becomes 


-n+l/n  r 

z  =  (h)  k  ^  /  rjdT  +  C 


C  +  ^  k  Z^idT 


n/n+1 


[Bl4] 
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Boundary  conditions 


^  At  T  =  1  j  I  =  1  ;  T]  =  1. 

U 


f 


> 

ii. 

M 

y 


Hence 


Hence  [Bl4]  becomes 


I  = 


^  n+l  -  ,  . 

1  +  -  k  /  TjdT 

n  ' 


n/n+1 


[B15] 
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If  n  =  2 


1  +  I  k 


T 

J  T)ClT 


When  Weibull  type  distributions  are  used. 


T)  =  1  -  exp  ( -T  ) 


[Bl6] 


At  T  =  1 


n  1  e-1 

T]!  =  1  -  —  =  - 

e  e 


ri  =  [  1-exp  (-t“)] 


dn  _  e-1  r  a-1 


dx  e 


[  ax  exp ( -X  )  ] 


At  X  =  1 


Hence 


dn  _  e 

dx  ^  ~  e-1 
x=l 


[  1-exp (-x“)] 


[B17l 


42 


From  Equation  [B17]^  I  is  calculated  as  a  function  of  t  for 
various  values  of  a  as  shown  In  Figure  19.  Figures  31  and  32 
show  the  same  relationships  when  n  =  3  and  4  respectively.  When 
Gumbel  type  distributions  are  used 

ri  =  exp[-|3  exp(-T)] 

r\i  =  exp[-  |] 

-  ^  exp[-g  exp(-T)] 

exp[  -jS/e] 


k  = 


_  ^ 


dr 


,  ■  expf-g/e]  f 

T=1 


k  =  iS/e 


Hence 


I  = 


1  +  I  I  / 


2/3 


[B18] 


Figure  33  shows  the  variation  of  I  with  t  for  three  values  of  j3. 
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Cumulative  Depth  of  Erosion 


IS-  .  . 

f  =  6  ^  e  .  ^  o  ,  fdr’  1 

I  dt  ■’  max  e  Idt) 

max  max  \  'max 


Boundary  conditions^ 


Then 


Hence 


tt  =  ti,  r  =  rij  r  =  r/rij  t  =  t/t. 


e  ri  dr 

I  ti  dT 

max 


S,r. 


dr  =  I  dx  +  C 


S,rx 


r  =  /  I  dT  +  C 


[BI9] 


Boundary  conditions 


t  =  0;  r  =  Oj  T  =  0;  1  =  0;  r  =  0 


Hence 


C  =  0. 
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Again  at 


t  -  ti;,  T  =  1,  f  =  1 


I  dT 


Equation  [BI9]  becomes 


I  dT 


I  dT 


[B20] 


The  variation  of  r  with  t  for  various  values  of  a  when  n  =  2  is 
shown  in  Figure  3^- 

Relationship  Between  the  Maximum  Intensity  of  Erosion  and 
The  Intensity  of  Collapse 

Prom  Equation  [Bll] 


-(n+1 )/n 


A  I  S 
c ,  e 


___  m.  ^ 

n-l/n  ti  dT 

'  T=1 
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Using  Welbull  type  distributions  and  the  value  n  =  2 


3/2 


max 


e-1 


a 


(A^I  )^S 
c  e 


ti(e-l) 


3/2 


max 


a 


1  S  A 
c  e 


2[e(e-l)] 


i  ti 


Relationship  Between  Erosion  Strength  and  Maximum  Rate  of 


From  Equation  [B2l] 


3/2  3/2  1/2  A  S 

(r  )  S  =  M  I 

max'  e  c 


where 


M  = 


a 


[e(e-l)] 


_  ^ 

max  dt 


at  t  =  tl 


max 


Hence 


A^I  if 


S  = 


max 


[B21] 


Erosion 


[B22] 
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TABLE  2 


Test  Results  and  Weibull  Distribution 
Function  for  Fatigue  Failures 


Material:  1100-0  Al. 

Sample  Size:  n  =  29  specimens 

Temperature:  Room  temperature 

water  bath 

Frequency:  l4.2  kcs 

Stress  Level:  5OOO  psi 

i 

i 

Mean  Rank 

Order 

No.  of  Cycles 

1  n  \ 

q 

to  Failure,  N 

100  ^  =  F(N)  X  100 

1 

2.14  X  10^ 

3.33 

2 

2.29  X  10^ 

6.66 

3 

4.00  X  10^ 

9.99 

4.27  X  10^ 

13.33 

5 

5.30  X  10^ 

16. 66 

6 

5.40  X  10^ 

19.99 

7 

9.25  X  10^ 

23.32 

8 

1.17  X  10^ 

26.65 

9 

1.93  X  10° 

29.99 

10 

4.49  X  10^ 

33.33 

11 

6.00  X  10^ 

36.66 

12 

8.64  X  10^ 

39.99 

13 

1.09  X  lo"^ 

43.33 

14 

1.92  X  10^ 

46.66 

15 

1.98  X  lo"^ 

49.99 

l6 

2.32  X  10^ 

53.33 

17 

2.56  X  lo"^ 

56. 66 

18 

2.56  X  10^ 

59.99 
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TABLE  2  (Continued) 


Order 

q 

No  of  Cycles 
to  Failure,  N 

19 

2.73  X  lo"^ 

20 

3.58  X  10^ 

21 

3.84  X  lo”^ 

22 

3.95  X  10^ 

23 

4.72  X  lo”^ 

24 

5.55  X  lo”^ 

25 

8.17  X  lo'^ 

26 

8. 99  X  lo"^ 

27 

1.075  X  10® 

28 

1.20  X  10^ 

29 

1.66  X  10^ 

Mean  Rank 

100  =  P(N)  X  100 

63.33 

66.66 

69.99 

73.33 

76.66 

79.99 

'•'  83 . 33 

86.66 

89.99 

93.33 

96.66 
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TABLE  3 


Test  Results  and  Welbull  Distribution  Function 
For  Fatigue  Failures 


Material:  3l6  . 

Sample  Size:  n  = 

Temperature:  Room 
Wate: 

S.  S,  Frequency:  14.2  kcs 

20  specimens  Stress  Level:  30^000  psi 

temperature 
r  Bath 

Order 

No.  of  Cycles 

Mean  Rank 

q 

to  Failure^  N 

100  =F(N)xlOO 

n+1  '  ' 

1 

1.04  X  10^ 

^.75 

2 

1.56  X  10^ 

9.50 

3 

1.85  X  10^ 

14.3 

4 

2.56  X  10^ 

19.0 

5 

3.37  X  10° 

23.8 

6 

4.70  X  10^ 

28.6 

7 

4.81  X  10^ 

33.3 

8 

5.53  X  10^ 

38.1 

9 

6.16  X  10^ 

43.0 

10 

1.37  X  lo"^ 

^7.5 

11 

1.39  X  lo"^ 

52.5 

12 

1.60  X  lo'^ 

57.0 

13 

5.97  X  lo"^ 

62.0 

14 

6.40  X  lo'^ 

66.6 

15 

1.02  X  10® 

71.5 

l6 

1.37  X  10® 

76.0 

17 

3.2  X  10® 

81.0 

18 

4.6  X  10® 

86.0 

19 

6.0  X  10® 

90.5 

20 

7.96  X  10® 

95.0 
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TABLE  4 


Experimental  Data  Showing  the  Effect  of  Impact 
Velocity  on  the  Rate  of  Volume  Loss 


316  Stainless  Steel:  Ui^i  =  I50  fps 

Rel.Vel. 

No. 

Velocity 

Peak  Rate 

(u-u  ) 

fps 

cc/hr 

min. 

fps^ 

1 

520 

20  X  10“3 

20 

370 

2 

436 

7  X  10“^ 

55 

286 

-4 

3 

350 

19  X  10 

120 

200 

4 

292 

8  X  10"^ 

300 

142 

1100-0 

Aluminum:  =  50  fps 

-4 

1 

200 

20  X  10 

25 

150 

2 

166 

-4 

11  X  10 

35 

116 

3 

182 

8  X  lo"^ 

12 

132 

4 

224 

20  X  10”^ 

3 

174 

5 

274 

31  X  10~^ 

2 

224 

54 


ROTATING  DISC 
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^  SPECIMEN  MOVING 


FIGURE  2  -  TYPICAL  WATER  COLUMN  AND  SPRAY  PATTERN  DURING 
SPECIMEN  PROGRESSION 
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FREQUENCY 

OSCILLOSCOPE  COUNTER 


FIGURE  3  -  BLOCK  DIAGRAM  OF  THE  MAGNETOSTRICTION  APPARATUS  USED 
FOR  HIGH  FREQUENCY  FATIGUE  TESTS 
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cb\w 


NOTE; 


1 .  GROOVE  (DETAIL  “  A" )  IS  TO  BE  SMOOTH,  FREE  FROM  CHATTER 
TOOL  MARKS,  GROOVES  OR  OTHER  DISCONTINUITIES.  THE 
DIMENSIONS  OF  THE  GROOVE  MUST  BE  IDENTICAL  FOR  ALL 
SPECIMENS  IN  A  LOT  +0.001  AS  MEASURED  WITH  AN  OPTICAL 
COMPARATOR . 

2.  FINISH  IN  GROOVE  ^  OR  BETTER. 

3.  DIMENSION  '•  A"  =  HALF  WAVE  LENGTH  OF  SOUND. 


FIGURE  4  -  NOTCHED  HIGH  FREQUENCY  FATIGUE  SPECIMEN 
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1  2 

—  AND  ARE  RELATIVE  LENGTHS 
A  A 


CURVE  OBTAINED  FROM  NEPPIRAS  THEORY 


X 


FIGURE  5  -  BASIC  APPROACH  FOR  THE  DESIGN  OF  DUMB-BELL 
SHAPED  FATIGUE  SPECIMEN 
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DECIMALS  ±0.001 

FRACTIONS  ±0.010 


NOTE:  0.625  r  REDUCED  AREA  TO  BE  GROUND 
AND  POLISHED  TO  16  RMS  MINIMUM 


FINISH  #  125  UNLESS  NOTED 


ALL  DIMENSIONS  ARE  IN  INCHES 


TYPICAL  ROOM  TEMPERATURE  DIMENSIONS 


MATERIAL 

DIMENSION  A 

DIMENSION  B 

316  STAINLESS  STEEL 

1.812 

4.375 

1100-0  ALUMINUM 

1  .855 

4.460 

FIGURE  6  -  DUMB-BELL  SHAPED  HIGH  FREQUENCY  FATIGUE  SPECIMEN 


61 


,01  X 


o 


CO  CM 


■Ni/'Ni  -  Nivais  aaansvaw 


X 


I 

O 


X 


o 


62 


FIGURE  7  -  COMPARISON  OF  THEORETICAL  STRAIN  AND  MEASURED  STRAIN  FOR 
DUMB-BELL  SHAPED  FATIGUE  SPECIMENS 


FIGURE  8  -  316  STAINLESS  STEEL  SPECIMEN  RUN  20  HOURS  AT  150  FT/SEC  FOR 
8.3  X  106  IMPACTS  TO  THRESHOLD 
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PREVIOUSLY  REPORTED 
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D3S/ld  -  AliD013A  N3WID3dS 
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IMPACTS  TO  INCEPTION 

FIGURE  10  -  RELATIONSHIP  BETWEEN  IMPACT  SPEED  AND  THE  NUMBER  OF  IMPACTS  UNTIL  INITIAL 
PLASTIC  DENTS  ARE  OBTAINED  ON  316  STAINLESS  STEEL  AND  ON  1100-0  ALUMINUM 
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FIGURE  11  -  1100-0  ALUMINUM  FATIGUE  DATA  AT  HIGH  FREQUENCY  -  RELATIONSHIP 
BETWEEN  STRESS  AMPLITUDE  AND  NUMBER  OF  CYCLES  TO  FRACTURE 


INTENSITY  OF  EROSION 
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304-L  STAINLESS  STEEL 
.0, 


MATERIAL: 

LIQUID:  WATER  AT  80°F 

FREQUENCY:  16  KCS 
SPECIMEN 

DIAMETER:  0.625  INCH 

DOUBLE 

AMPLITUDE:  A  2.50  X  lO”^  INCHES 
O  1.85X  10”^  INCHES 
O  1 .25  X  10'^  INCHES 
□  0.65  X  10*^  INCHES  H 


0  200  400  600  800  1000  1200  1400  1600  1800  2000  2200 

TIME -MINUTES 

FIGURE  13  -  EFFECT  OF  TIME  ON  CAVITATION  DAMAGE  RATE 
FOR  VARIOUS  AMPLITUDES  (REFERENCE  9) 


FIGURE  14  -  EFFECT  OF  TIME  ON  INTENSITY  OF  EROSION 
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MATERIAL  316  Stainless  Steel 
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FIGURE  15  -  RATE  OF  VOLUME  LOSS  AS  A  FUNCTION  OF  TIME  AT  THREE  VELOCITIES 
FOR  316  STAINLESS  STEEL 


SPECIMEN  SPECIMEN 
VELOCITY  FT/SEC  NO. 
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MATERIAL  1100  -  0  Aluminum 

SPECIMEN  cDcr'iJiiCKi 

DENSITY  2.71  g/cc  VELOCITY  SYMBOL 
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FIGURE  17  -  RATE  OF  VOLUME  LOSS  AS  A  FUNCTION  OF  TIME  AT  THREE  VELOCITIES 
FOR  1100-0  ALUMINUM 


VELOCITY  SPECIMEN 
FT/SEC 


FIGURE  19  -  THEORETICAL  PREDICTION  OF  THE  EFFECT  OF  TIME 
ON  INTENSITY  OF  EROSION  WHEN  n  =  2 


IMPACT  VELOCITY 
FT/SEC 


FIGURE  20  -  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  THE  EROSION  THEORY 
FOR  1100-0  ALUMINUM 


IMPACT 

VELOCITY  FT/SEC 
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FIGURE  21  -  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  THE  EROSION  THEORY 
FOR  316  STAINLESS  STEEL 


PER  CENT  FAILED  -  F  (N)  x  100 


MATERIAL 
SAMPLE  SIZE 
FREQUENCY  OF  TEST 


1100-0  ALUMINUM 
29 

14.2  kcs 


STRESS  LEVEL 
TEMPERATURE 
ENVIRONMENT 


5000  psi  (Axial  Push-Pull) 
ROOM  TEMPERATURE 
WATER  COOLED 


NUMBER  OF  CYCLES  TO  FAILURE  -  N 


FIGURE  22  -  WEIBULL  DISTRIBUTION  FOR  HIGH  FREQUENCY  FATIGUE  OF 
1100-0  ALUMINUM  AT  ROOM  TEMPERATURE 
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PER  CENT  FAILED  -  F  (N  )  x  100 


MATERIAL 
SAMPLE  SIZE 
FREQUENCY  OF  TEST 


316  STAINLESS  STEEL 
20 

14.2  kcs 


STRESS  LEVEL 
TEMPERATURE 
ENVIRONMENT 


30,000ps1  (Axial  Push-Pull) 
ROOM  TEMPERATURE 
WATER  COOLED 


NUMBER  OF  CYCLES  TO  FAILURE  -  N 

FIGURE  23  -  WEIBULL  DISTRIBUTION  FOR  HIGH  FREQUENCY  FATIGUE  OF 
316  STAINLESS  STEEL  AT  ROOM  TEMPERATURE 
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PEAK  RATE  OF  VOLUME  LOSS  -  cc 


8  n  1/1 6  I NCH  DIAMETER  WATER  JET' 


11 00-0  ALUMINUM 


316  STAINLESS  STEEL 


40  60  80  100 


400  600  800  1000  2C 


IMPACT  VELOCITY  -  fps 


FIGURE  24  -  RELATIONSHIP  BETWEEN  PEAK  RATE  OF  VOLUME  LOSS  AND 
THE  IMPACT  VELOCITY 


TIME  CORRESPONDING  TO  PEAK  RATE  OF  EROSION,  f,  -  MINUTES 


IMPACT  VELOCITY  -  fps 


FIGURE  25  -  RELATIONSHIP  BETWEEN  THE  TIME  AT  WHICH  PEAK  RATE  OF 
VOLUME  LOSS  IS  OBSERVED  AND  THE  IMPACT  VELOCITY 
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PEAK  RATE  OF  VOLUME  LOSS  -  ccAr 


RELATIVE  IMPACT  VELOCITY  (U-U^)  -  fps 


FIGURE  26  -  RELATIONSHIP  BETWEEN  THE  RELATIVE  IMPACT  VELOCITY  AND 
THE  RATE  OF  VOLUME  LOSS 
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1  mm 


I 


nOO-O  ALUMINUM 


-  1  mm 


316  STAINLESS  STEEL 


FIGURE  27  -  MICROPHOTOGRAPHS  OF  THE  TEST  MATERIALS 
SHOWING  STRUCTURE 
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LIQUID  IMPACT  EROSION  OF  1100-0  ALUMINUM 
IMPINGEMENT  VELOCITY  450  fps  FOR  3  MINUTES 
(MAGNIFICATION:  88X  ) 


CAVITATION  EROSION  OF  1100-0  ALUMINUM  USING  THE 
MAGNETOSTRICTION  DEVICE  DRIVING  FREQUENCY 

14.2  kcs  DOUBLE  AMPLITUDE  2.0  x  lO"^  INCHES  FOR 
1  MINUTE  (  MAGNIFICATION:  74X  ) 


FIGURE  28  -  COMPARISON  OF  INDENTATIONS  PRODUCED  BY  MULTIPLE 
LIQUID  IMPACT  EROSION  AND  BY  CAVITATION  EROSION  - 
SCANNING  ELECTRON  MICROPHOTOGRAPHS 
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FIGURE  29  -  ROTATING  DISC  CAPABLE  OF  HOLDING  SIX  SPECIMENS 
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RELATIVE  INTENSITY  -  I  RELATIVE  INTENSITY 


RELATIVE  TIME  -  T 

FIGURE  33  -  THEORETICAL  PREDICTIONS  OF  THE  EROSION  INTENSITY  WHEN 
GUMBEL  TYPE  DISTRIBUTIONS  ARE  USED 


0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0 

RELATIVE  TIME  -  T 


FIGURE  34  -  THEORETICAL  PREDICTIONS  OF  THE  CUMULATIVE 
DEPTH  OF  EROSION 
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NASA-Langley,  1969 
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